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Ferromagnetism in Mn-doped GaAs, the prototypical dilute magnetic semiconductor, has so far 
been attributed to hole mediated RKKY-type interactions. First-principles calculations reveal a 
strong direction dependence of the ferromagnetic (FM) stabilization energy of two magnetic ions, 
a dependence that cannot be explained within RKKY. In the limit of host-like hole (engineered 
here by an GGA-f U approach with large U) where the RKKY model is applicable, we find that 
the exchange energies are strongly reduced, suggesting that this limit cannot explain the observed 
ferromagnetism. The dominant contribution stabilizing the FM state is found to be maximal for 
<110>-oriented pairs and minimal for <100> oriented pairs, providing an alternate explanation 
for magnetism in such materials in terms of energy lowering due to p-d hopping interactions, and 
offering a new design degree of freedom to enhance FM. 

PACS numbers: PACS number: 75.50Pp, 75.30.Hx 



The discovery of ferromagnetism in Mn doped GaAs 
has spurred considerable attention in this important 
class of materials. Experimentally it is known that the 
introduction of Mn in GaAs gives rise to an acceptor . 
The hole produced by the acceptor is believed to inter- 
act with the localized orbitals of the TM impurity and 
mediate ferromagnetism. The question is what type of 
mechanism explains the FM resulting from the interac- 
tion between the hole and the magnetic ion. In a model 
Hamiltonian approach ^ |^ one selects a priori a fa- 
vored mechanism and works out its physical consequences 
and manifestations. For example, in the limit where the 
magnetic electrons can be treated as a localized entity, 
and the quantum oscillations of the electron spin polar- 
ization around the localized impurity can be neglected, 
the exchange interaction between the TM impurity and 
the hole can be RKKY-like. It has been argued IJ] that 
this limit is indeed reached for TM impurities in semicon- 
ductors. A consequence is that the exchange interaction 
has either a vanishing dependence on the direction of the 
vector joining the Mn ions in GaAs (if a spherical Fermi 
surface is assumed), or a weak one ;5] if the true non- 
spherical fermi surface of the host is considered. 

As an alternative one can use ab-initio total energy cal- 
culations for magnetic ions in a host crystal to distill a 
mechanism a posteriori. We consider TM (V-Fe) pairs in 
GaAs, at various separations and calculate the exchange 
interaction strength, Jy (R). For all cases Jy (R) are 
found to exhibit a strong dependence on the specific lat- 
tice orientation of the TM pairs, in sharp contrast to the 
simplest realization of the RKKY model with a system- 
independent spherical Fermi surface. To test if an ex- 
tended RKKY model does better, we have calculated the 
anisotropic Jrkky(R) Jj, taking the Fermi surface of 
hole doped GaAs explicitly into account. We find that 



Jrkky(^) is qualitatively different from Jy(R) deter- 
mined from ab-initio calculations, thereby establishing 
that the magnetic interactions in these systems cannot 
be described even within a realistic RKKY-type model. 
The ab-initio results are subject to specific uncertainties 
in the energy position of the d levels [sl • To see if this can 
affect our conclusion we use a simplified self-interaction 
correction scheme in the form of GGA-I-U 9]. We tune 
U so as to fit the incorrect GGA value of the energy po- 
sition of the primarily Mn d states in the valence band of 
GaAs {E.,j - 2.6 eV) to experimental photoemission (Ey 
-4 eV) |l(|. The strong non-RKKY anisotropy is still 
present for U ^ 3-4 eV, proving that the GGA error is 
qualitatively inconsequential. Finally we show that this 
directional dependence can be explained within a model 
of ferromagnetism arising from energy gain coming from 
p-d hopping interactions 11] 

We have carried out first-principle electronic structure 
calculations using density functional theory, within the 
pseudopotential plane- wave total ener gy m ethod [l^ , us- 
ing ultra-soft pseudopotentials (USP) [l^l and projected 
augmented wave (PAW) [13 potentials as implemented 
in VASP code [ij. The equilibrium lattice constant of 
the TM containing GaAs supercells was fixed at the value 
obtained for zincblende GaAs (5.728 A) using the PW91 
GGA exchange functional 0|, but the atomic positions 
were allowed to relax. The basis sets had a cutoff energy 
for plane waves equal to 13.3 Ry. K-point grids of 4x4x4 
including F and 2x1x1 were used for the 64 and 256 atom 
calculations. GGA-I-U calculations were performed with 
a J7 on Mn, keeping the intra-atomic exchange interac- 
tion fixed at values used earlier [l^ . while U was varied. 

In order to understand the role played by the hole, we 
consider the cases of V and Fe in GaAs, both of which do 
not introduce holes into the system. Fig. 1(a), (b) show 
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the TM d projected partial density of states (PDOS) re- 
solved into t2 and e symmetries for up (+) and down 
(-) spin channels. In each spin channel we have a pair of 
states (bonding and antibonding) with t2 symmetry. The 
magnetic ground state that would be favored can be read- 
ily understood with a schematic two level model shown 
in Figs. 2(a) and (b). The unperturbed exchange-split 
3d levels on the isolated atoms TMl and TM2 are shown 
on the left and right side of Fig. 2(a) and Fig. 2(b) for 
FM and AFM (antiferromagnetic) arrangement of TM 
spins, respectively. The up and down spin states on the 
TM atoms interact via spin-conserving hopping interac- 
tions of strength v and form a set of bonding-antibonding 
states for each spin channel, as shown in the central part 
of each panel. In a FM arrangement (Fig. 2(a)), both 
bonding and antibonding levels of one spin channel are 
completely filled, so to a first order, there is no gain in 
energy in this magnetic coupling. For the AFM arrange- 
ment (Fig. 2(b)), however, the bonding states with t2 
symmetry are completely filled for both spin channels, 
while the antibonding states are empty. Consequently, 
the resulting AFM energy gain is ~ v"^/!, where / is the 
energy separation of the same spin levels on TMl and 
TM2. Hence, the AFM arrangement of the TM spins is 
favored in the absence of a hole. The expectations of the 
simple model of Fig. 2 are verified by the results from 
our ab-initio calculations (Fig 3(a) and (b)). The AFM 
configuration is favored at all separations, with the excep- 
tion of V at first neighbor. Interestingly the largest AFM 
stabilization energy is only 31 meV for V, while it is 298 
meV for Fe. This difference can be understood in terms 
of the hopping interaction strength, v, entering the u^// 
stabilization of the AFM states. When the highest occu- 
pied states have t2 symmetry as in GaAs:Fe (Fig. 1(b)), 
the relevant hopping matrix element is between the Fe 
t2 states. These are much larger than those between e 
states as in GaAs:V (Fig. 1(a)) because e{t2g) orbitals 
point in-between (towards) the nearest-neighbors. 

Turning next to GaAs:Mn and GaAs:Cr, it is evident 
from the PDOS (Figs. 1 (c) and (d)) that both these 
impurities introduce holes in the system. In the presence 
of partially occupied orbitals, the simple model of Fig. 
2 predicts ferromagnetism as the energy gain for a FM 
arrangement is large because the interacting levels are 
degenerate in the case of FM arrangement, while these 
are separated by a large energy in the AFM case. The 
expectations of the simple model are verified by our ab- 
initio calculations (Figs. 3(c) and (d)). Ferromagnetism 
is favored at all separations for Cr and Mn pairs. 

Focusing on Mn-doped GaAs, we extract Jij from 
FipM-^AFM of Fig. 3(d) for different orientations of Mn 
atoms in the 64-atom cell, as well as for the 256 atom 
cell. The significant feature of Jy shown in Figs. 4(a) 
and (b) is the pronounced domination of orientation over 
distance dependence. In Fig. 4(a) the three pairs ori- 
ented along the <110> direction (connected by a dotted 



line) show a monotonic decay with R, while remaining 
higher in strength compared to the pairs oriented along 
other directions {e.g. <100> direction, connected by a 
dashed line), even when such pairs have a smaller sepa- 
ration. This is further established by our results for two 
Mn atoms at the same distance, but oriented in differ- 
ent directions, namely <110> and <411>. One Mn is 
placed at the origin and the other either at (1.5a 1.5a 0) 
for <110> or at (2a 0.5a 0.5a) for <411>. The calcu- 
lated Jij's for these two pairs at the same separation 
are vastly different (Fig. 4(a)). Such an observation 
is obviously incompatible with the usual RKKY model 
based on an isotropic Fermi surface. It is however, pos- 
sible that such orientation dependencies arise from the 
non-spherical Fermi surface of the specific system. We 
have calculated the orientation dependent exchange in- 
teraction strengths, Jrkky based on the RKKY model 
including the realistic band structure effects such as the 
non-spherical Fermi surface of the host GaAs. The 64 
atom supercell of GaAs with one hole was taken and 
the eigenvalues were computed over a grid of 6x6x6 k- 
points. The eigenvalues were interpolated over a finer 
grid of 10x10x10 and the generalized susceptibility x(q) 
was computed using the method of Ref . 01 . The Fourier 
transform of xil) was used to calculate Jrkky- We 
checked the stability of our calculation by increasing the 
number of k-points to 20x20x20. The changes were found 
to be less than 5%. This Jrkky is plotted for compari- 
son as an insert to Fig. 4(b). Evidently, the behaviors of 
Jij and Jrkky are qualitatively different; for example, 
the first principles calculated Jy is smallest along <100> 
and largest along <110> as seen in Fig. 4, whereas 
Jrkky is almost maximal for <100>. Obviously, any 
RKKY-type model in spite of extending it to account for 
real band structure effects is inadequate to describe the 
underlying magnetic interactions of these systems. 

The above mentioned failure of RKKY model is in 
fact easy to understand, as GaAs:Mn clearly violates the 
fundamental assumptions needed for the validity of the 
RKKY model. The RKKY theory involves a perturba- 
tivc treatment in which the exchange splitting (Eexch) 
of the host band is small in comparison with the Fermi 
energy (Ep), Eexch << Ep. However, the DMS's, in 
particular Mn doped GaAs, are half-metallic ferromag- 
nets, with complete spin-polarization which arises from 
Eexch being larger than Ep. Thus, a perturbation in 
Eexch/ Ep is bound to fail, making the inapplicability of 
RKKY mechanism obvious for these systems. Another 
interesting consequence of the half-metallicity is the com- 
plete supression of spin flip scattering between up and 
down spin states of the conduction electrons essential in 
the RKKY exchange coupling, thereby distinguishing the 
present system from those dominated by RKKY interac- 
tions. It should be noted that total Jrkky is a prod- 
uct of two terms. The first term is proportional to the 
square of the strength of the spin-coupling between the 
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local (Mn) moment and the conduction electrons explic- 
itly accounted for in the Kondo-lattce Hamiltonian; the 
second term includes all the band structure information 
concerning the host lattice. All RKKY-type approaches 
assume the first term to be a constant, representing the 
strength of the spin-coupling between the local moments; 
thus, all the dependencies on the distance and orientation 
within RKKY approach arise exclusively from the second 
term. We have already shown that the R dependence of 
Jrkky in the inset to Fig. 4(b) is entirely inadequate 
to describe the Jij (R) observed. Next we point out that 
the R dependence of is in fact controlled almost en- 
tirely by the distance and the orientation dependencies 
of the spin-coupling in the Kondo-lattice model, which 
itself arises from the anisotropic hopping for example in 
a Periodic Anderson Hamiltonian. 

A single Mn in GaAs introduces fully occupied t+, e+ 
states inside the valence band, and partially occupied i+ 
state at Ep made of TM d and anion p orbitals. These 
partially occupied levels are represented in the left and 
right panels of Figs. 2(c) and (d). They interact via 
hopping and lower the total energy of the FM arrange- 
ment. The dependence of the exchange integral on lat- 
tice orientation comes from the dependence of the hop- 
ping matrix element entering the FM energy stabiliza- 
tion. This is different from any dependencies within the 
RKKY mechanism that arise from non-spherical Fermi 
surface 19]. The mechanism discussed here based on p-d 
hopping is not unique to dilute magnetic semiconductors, 
but is common to a wide class of materials. It was first 
introduced to explain the robust ferromagnetic state of 
Sr2FeMo06 0|. In the present work, we have pointed 
out another novel aspect of this mechanism in terms of 
its specific and characteristic orientation dependence. 

It is interesting to examine whether the orientation 
dependence changes with the localization of the hole- 
carrying orbital. We achieve this using the GGA+U 
approach (Si] with a finite that pushes the bonding 
levels at Ev-2.6 eV (Fig. 1(d)) deeper in the GaAs va- 
lence band, making them more Mn-localized, while the 
hole-carrying state at Ep becomes more host-like and 
delocalized. Figs. 5(a) and (b) show the Mn d PDOS 
with t2 symmetry for [/ = 0, 6, 10 and 15 eV. As is ev- 
ident from the inset of Fig. 5(a), the introduction of U 
pushes the location of the Mn feature from Ey-2.6 eV at 
U =0 to Ey-5, Ey-7 and Ey-9.5 eV for U =6, 10 and 15 
eV, respectively. Agreement with the photoemission de- 
termined position [ill of £"„ - 4 eV requires a f7 of around 
3-4 eV. Most features of the U=0 calculations are pre- 
served at this value of U, including the strong anisotropy 
in Jij (see Fig. 4(b)). Thus the GGA error does not 
affect our results much. 

We can use GGA-I-U to simulate the condtions under 
which RKKY is supposed to work: The amplitude of 
the Mn d PDOS of the anti-bonding i+ states at Ep de- 
creases as U increases (Fig. 5 (a),(b)). This decrease in 



Mn content is clearer from the hole wavefunction squared 
plotted in the <110> plane for U =0 and 10 eV in Figs 
5 (c) and (d): At U =0, a considerable portion of the 
hole wavefunction at Ep is localized on Mn and its near- 
est neighbor As atoms, while at C/ =10 eV, the states at 
Ep become more delocalized, host-like as in the case for 
GaAs:Zn. At this limit (;7=10-15 eV) of " host-like-hole" 
the conventional RKKY approach is supposed to be valid. 
Our calculations show that at this limit the FM stabiliza- 
tion J is already quite small, and the Jy 's become more 
short-ranged with only nearest-neighbor pairs contribut- 
ing (Fig. 4(b). Thus, the observed FM is unexplained by 
a model simulating "host-like-hole" RKKY conditions. 

In summary, we have examined the microscopic mech- 
anism giving rise to ferromagnetism in 3d impurities in 
GaAs. A strong deviation is found from current carrier- 
mediated ferromagnetism based models 0, , which we 
find arc not appropriate even when the hole is more host- 
like. The dominant contribution to FM stabilization is 
found to be from p-d hopping. 
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FIG. 1: (Color online) The broadened up(+), down(-) spin 
TM d PDOS in spheres of radius 1.2 A with t2, e symmetry 
for different TMs. 



FIG. 2: (Color online) Schematic energy levels for two inter- 
acting TM with their spins FM [(a),(c)] and AFM [(b), (d)] 
aligned and highest occupied level fully [(a), (b)], partially 
[(c),(d)] filled. 



FIG. 3: Distance/Orientation dependence of Efm/^^afm for 
two (a) V, (b) Fe, (c) Cr , (d) Mn in 64 atom GaAs cell 
using USP potentials (using PAW in (d) in parentheses). The 
upper a;-axis gives the direction of the vector joining the two 
TM atoms. 



FIG. 4: The distance/orientation dependence of Jij for Mn 
pairs in (a) 256, (b) 64 atom GaAs cell using PAW potentials. 
The expected dependence of Jrkky for a hole in GaAs is 
given in the insert. 



FIG. 5: ( Color) The up [(a) and inset] and down (b) spin 
Mn t2 PDOS for U=0(black line), (7=6(red line), C/=10(green 
line) and f/=15 (blue line) eV. Hole wavefunction squared in 
the <110> plane are shown for ?7=0, and ?7=10 in parts (c) 
and (d) respectively. 
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